Anderson localization in the quintic nonlinear Schrodinger equation 
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In the present paper we consider the quintic defocusing nonlinear Schrodinger equation in presence of a 
disordered random potential and we analyze the effects of the quintic nonlinearity on the Anderson localization 
of the solution. The main result shows that Anderson localization requires a cutoff on the value of the parameter 
which controls the quintic nonlinearity, with the cutoff depending on the amplitude of the random potential. 
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Introduction - The Anderson localization (AL) phe- 
nomenon is the suppression of transport due to a destructive 
interference of the many paths associated with coherent mul- 
tiple scattering from the modulations of a disordered potential 
[1]. This effect has been experimentally observed for light in 
diffusive media [2, 3], photonic crystals [4, 5], optical fiber 
arrays [6], microwaves [7], sound waves [8], and others. Re- 
cently, AL was also observed in noninteracting Bose-Einstein 
condensate (BEC) [9, 10]. For instance, in Ref. ['-)] the au- 
thors used a 87 Rb BEC in a one-dimentional (ID) waveguide 
in the presence of a controlled disordered potential generated 
by a laser speckle for observation of an exponential tail of the 
spatial density distribution, which is a signature of the AL; 
the AL was also observed in [10] for a 39 K BEC in presence 
of a ID quasiperiodic bichromatic optical lattice. Motivated 
by these experimental investigations, many theoretical studies 
has been done in the last years considering random potentials 
[1 1-15] and bichromatic optical lattice [16]. 

In BECs, the search for localized structures such as solitons 
and breathers has also attracted attention of many researchers. 
As example, we have the experimental observation of dark 
solitons [17], formed in BEC with repulsive 87 Rb atoms, and 
bright solitons [18], generated in attractive 7 Li atoms. From 
the theoretical point of view, the control of such solutions 
can be facilitated through the search for analytical solutions 
of the ID Gross-Pitaevskii equation (GPE). In this sense, re- 
cently, analytical solitonic solutions to the more general case, 
employing space- and time-dependent coefficients, was con- 
sidered for the cubic [19], the cubic-quintic [ 10], the quintic 
[21], and also the GPE in higher dimensions [22]. Analytical 
breather solutions has been found in Ref. [23]. 

A current challenging problem in BECs is to understand 
how the nonlinear effects may affect the localization which 
appears from the presence of random potential. In this 
sense, the destruction of AL by a weak nonlinearity was re- 
cently studied considering a one-dimensional discrete nonlin- 
ear Schrodinger lattice with disorder [14]. In Ref. [14], the 
authors demonstrated numerically that above a certain criti- 
cal strength of cubic nonlinearity the AL is destroyed and an 
unlimited subdiffusive spreading of the field along the lattice 
occurs. Here, we consider a quintic nonlinear Schrodinger 
equation (QNLSE) in presence of a weak random amplitude 
potential. Our goal is to study the effects of the quintic repul- 
sive (defocussing) nonlinearity on the AL of the system. In 



the BEC case, the quintic nonlinearity is related to the three- 
body scattering [ ], or to a mere expansion of the cubic term 
due to a reduction of dimensionality [25], while in nonlinear 
crystals it is due to the quintic order nonlinear electric suscep- 
tibility [26]. Recently, we considered the unidimensional 
reduction of a three-dimensional BEC with two- and three- 
body interactions [27], and there it was shown that the quintic 
GPE is the effective equation that describes the profile of the 
BEC when it presents a weak coupling regime, with the cubic 
nonlinearity being much smaller than the quintic one. 

Theoretical model - We deal with a ID quintic NLSE writ- 
ten in its dimensionless form [27] 



iipt = - \i>2 



(i) 



where ip = 4>{x, t) is the wave-function describing the collec- 
tive state of the atoms in a BEC or the electric field propagat- 
ing in a nonlinear crystal, ip t = dip/dt and ip xx = d 2 i\)jdx 2 , 
V = V(x) is the potential, and 95 is the nonlinear constant 
which controls the quintic nonlinearity. 

In the present paper we work with random potentials that 
are generated by an algorithm of random number of the 
MAPLE 13 program called RAND. Here, we use such algo- 
rithm to obtain a pattern of potential that present a random 
behavior, given by 



V(x) = V £(x), 



(2) 



with £(2;) = An for a; G [-L + 2nL/S, -L + 2(n + l)L/S] 
and n = 0, 1, 5—1, where A n is a random number, 2L 
the dimension of the system in spatial x direction and S the 
number of disordered amplitude points. Also, Vo is an ampli- 
tude parameter of the random potential in (2). We stress that 
the Eq. (2) consists of a perturbation in the amplitude, in the 
spatial range [-L, L~\. 

Numerical results - The Eq. (1) is numerically solved using 
a split-step algorithm based on the Crank-Nicolson method. 
Here we use space and time steps given by 0.04 and 0.001, 
respectively. First, we employ an imaginary-time method 
for propagation of a Gaussian pulse, to get the profile of the 
ground state solution. Next we realize various analyzes for the 
conditions for occurrence of localized solutions in the system. 
The interval of the system used in the simulations is [-L, L] 
with L = 30. We also selected S = 300, that is a typical 
number of perturbed points for the considered range [15]. 
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FIG. 1 : (Color online) Random amplitude potential given by Eq. (2), 
for Vb = 1. 




FIG. 2: (Color online) Profile \ij>\ 2 of the QNLSE in the presence of a 
potential with random amplitude considering g& — in solid (black) 
line, g$ = 1 in dashed (red) line, and gs — 3 in dash-dotted (green) 
line. Also, we have used the values Vo — 5 and S = 300. The arrow 
(blue) shows the position in which the fragmentation starts. 



In Fig. 1 we display the potential with random amplitudes 
given by Eq. (2). In BECs the random potentials can be mod- 
eled by laser speckles [9, 28] that is formed when a laser beam 
crosses a diffusive plate, and focused onto the atomic sample, 
creating a random phase shift along the beam profile, which 
is converted into a random intensity distribution. Also, the 
disorder can be produced by quasiperiodic lattices [10] that 
is generated through combination of a deep primary lattice 
with a secondary one via the dipole potential of a laser beam 
arranged in standing wave configuration. This quasiperiodic 
lattices also appears in quasicrystals [29]. 

Next we use the simulations to verify the influence of the 
quintic term of the QNLSE on the localization. Fig. 2 shows 
the stationary state obtained via imaginary time evolution of 
the QNLSE for three distinct values of the nonlinearity. The 
arrow (blue in online version) shows the position of a second 
peak associated to the fragmentation of the solution. This re- 
sult will be clarified below. 

We have analyzed the average (x) and peak x p positions 
of the solutions as well as the peak height (max \ip\ 2 ) and its 
derivative (dmax \ip\ 2 ) versus g$ as a tentative of verify the 
influence of the nonlinear quintic term. The results are plotted 
in Fig. 3a and 3b, considering Vq = 5. We look the lack of co- 
incidence in this values as indicative of the fragmentation in a 
multi-peak solution. Also, in Fig. 4 we exhibit the root-mean- 
square width Ax of the profile versus the nonlinear coefficient 
35 for different values of the amplitude of the random poten- 
tial. Note that an abrupt change in Ax occurs for a given value 
of the nonlinearity. This result is more visible for large val- 
ues of the amplitude of the random potential. For example, 
considering the case of Vq =3 (boxes) the abrupt transition 
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FIG. 3: (Color online) Influence of the nonlinear quintic term on 
(a) the average (a;) (circles) and peak x p (boxes) positions and (b) 
the peak height (max \ip\ 2 , circles) and its derivative (dmax \ip\ 2 , 
boxes). A lack of coincidence in the values plotted in (a) indicates 
a fragmentation in a multi-peak solution. We have used the values 
V = 5 and S = 300. 
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FIG. 4: (Color online) Root-mean-square width of the profile (Ax) 
versus the nonlinear coefficient gs . We display the case with Vo = 1 
in circles (yellow), Vo — 3 in boxes (red), and Vb = 5 in diamonds 
(green). 



appears at 35 ~ 1.2, and we hope to obtain exponential tails 
for values of 35 < 1.2. 

The influence of the amplitude of the random potential on 
Ax and its derivative is displayed in Figs. 5a and 5b, con- 
sidering #5 = and 95 = 1, respectively. Note that for the 
case with 35 = 0, Ax is decreased when Vq is increased. This 
result is expected since for large values of the amplitude of 
the potential the state tends to become strongly localized [15]. 
On the other hand, in the presence of a quintic nonlinearity 
Ax suffers abrupt changes. Then, only for a 'large' amplitude 
Vb the state turns localized (Vb ^ 5). 

We have search regions with localized solutions present- 
ing pure exponential decaying tails [ ■■']. Due to the asym- 
metry of the solution, we employed exponential fit (\tp\ 2 = 
TVexp[— 2\x\/l], where TV is the amplitude) for the left and 
right tails. The localized state presents a pure Gaussian tail 
when the system engender a strong random amplitude poten- 
tial, since in this case I — > 0. However, the exponential tail 
is obtained and the limit of localization in a weak random po- 
tential is attained when I > Ax [9, 10]. In this sense, in Table 
I we display comparative values of Ax, left tail II, and right 
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FIG. 5: (Color online) Root-mean-square width of the profile (Ax, 
circles) and its derivative (Ax', boxes) versus the potential amplitude 
Vo. In (a) we display the case with g$ = while in (b) it is shown 
the case with g$ = 1. 
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TABLE I: Comparative values of the root mean square width Ax, 
left tail II, and right tail Ir versus the potential amplitude Vo of the 
potential with random amplitudes considering two values of nonlin- 
earity (35 = and g 5 — 1.0). 
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FIG. 6: (Color online) Example of a profile with AL solution, consid- 
ering Vb = 5 and gs = 1. In solid (green) line we exhibit the profile 
\ip\ 2 of the solution while in dash-dot (blue) line it is shown the ex- 
ponential fit (~ exp(— 2|x — x p \/l)). For comparison, we display a 
Gaussian fit in dashed (red) line. 
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FIG. 7: (Color online) Peak height (max IV'I 2 ) vs. the number of 
disturbed points S, considering Vo = 1 with = (circles, yellow) 
and gs = 1 (boxes, red). 



analyzed the effects of the quintic nonlinearity on the Ander- 
son localization of the solution. A cutoff value for the quintic 
term is obtained related with the amplitude of the random po- 
tential, as well as that to the nonlinear Schrodinger equation 
with cubic nonlinearity. The method is general and can be 
applied for other configurations of disordered potentials. 
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tail In versus Vq considering two specific values for the non- 
linearity (<?5 = and 95 = 1.0). In the same way, Fig. 6 
shows an example of a profile with AL solution (solid line), 
the Gaussian (dashed), and exponential tail (dash-dotted) fits, 
considering Vq = 5 and 95 = 1. 

Finally, in Fig. 7 we display a simulation of the peak height 
(max lip] 2 ) versus the number of disturbed points S, consider- 
ing Vq = 1 with 95 = (circles) and 95 = 1 (boxes). Note 
that in the region S < 200 the peak height oscillate and for 
S > 200 this value is 'stabilized'. So, we justify the use of 
the S = 300 disturbed points for the potential considered in 
the present simulations such that a weak change in this value 
does not interfere in the results. To this end, we have tested 
the accuracy of the numerical simulation by varying the pat- 
tern of the random amplitude potential in Eq. (2) and the total 
number of time and space steps. 

Conclusion - In conclusion, in the present paper we have 
consider the quintic defocusing nonlinear Schrodinger equa- 
tion in presence of a disordered random potential and we have 
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